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Quantifying experimentally created entanglement could in principle be accomplished by measuring 
the entire density matrix and calculating an entanglement measure of choice thereafter. Due to 
the tensor-structure of the Hilbert space, this approach becomes infeasible even for medium-size 
systems. Here we present methods to quantify the entanglement of arbitrarily large two-colorable 
graph states from simple measurements. The measurement data considered here is merely given by 
stabilizer measurements, thus leading to an exponential reduction in the number of measurements 
required. We provide analytical results for the robustness of entanglement and the relative entropy 
of entanglement. 
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Introduction Detecting [T] and quantifying entangle- 
ment [5] is one of the major tasks in quantum informa- 
tion science. Experimentally created entanglement can 
in principle be quantified by determining the quantum 
state via full tomography, and calculating an entangle- 
ment measure of choice for this state. Apart from ex- 
ceptions like the negativity, entanglement measures usu- 
ally involve optimization problems, which makes them 
hard to calculate. Another issue is the tensor-structure 
of the Hilbert space, which implies that the number of 
measurement settings grows exponentially with the num- 
ber of constituents involved in the system. Despite the 
recent developments in efficient tomography [3], the de- 
termination of the full quantum state appears to involve 
an unnecessary overhead given that only a single num- 
ber, the value of the entanglement measure, is required. 
For this reason, more sophisticated methods for the direct 
quantification of entanglement in many-body systems are 
required. 

Here we present direct and experimentally efficient 
methods to quantify entanglement of quantum many- 
body systems. We put the emphasis on two-colorable 
graph states, which represent a vast resource for ap- 
plications in quantum information science. They en- 
compass Greenberger-Horne-Zeilinger (GHZ) states [4], 
Calderbank-Shor-Steane (CSS) error correction code- 
word states, and cluster states [3]. Due to the impor- 
tance of graph states, a considerable experimental effort 
has been made to realize them using photons [6- ICS], cold 
atoms [11J . and proposals for trapped ions are pursued 
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We will show that the entanglement - according to 
a variety of entanglement measures [23 - of such two- 
colorable graph states can be estimated efficiently via 
measurements of the stabilizer operators only, thus re- 



ducing the experimental effort in measuring the state ex- 
ponentially. Furthermore our method of entanglement 
estimation is purely analytic, thus avoiding computation- 
ally costly post-processing of measurement data. 

Entanglement estimation Graph states of n qubits 
correspond to a graph G of n vertices, with n binary 
indices (ki, . . . ,k n ). We denote the Pauli matrices at 
the i-th qubit by Xi , , Z± . One can show that the 
2 n graph states |Gfc ll ....fc„) are the simultaneous eigen- 
states of the n mutually commuting operators: Ki := 
^i®Ngb(i)^j-i * = 1, where Ngb denotes the 

set of all neighbours of qubit i defined by the graph. 
Graph states satisfy the following eigenvalue equation: 
Ki\G kl ,...,k n ) = (— l) fei |C?fc 1 ,...,fc„). The n operators K l 
generate an abelian group S, called the stabilizer. An 
experimentally created graph state could in principle be 
verified by measuring the 2™ elements of the stabilizer. 
As mentioned in the introduction, full-state tomogra- 
phy is not an option for determining the properties of a 
quantum many-body system due to an exponentially fast 
growing measurement effort. We will see that merely the 
measurement results of the generators of the stabilizer 
suffice to attain highly useful bounds on entanglement 
measures. 

Let us suppose the goal of an experiment is the creation 
of a two-colorable graph state, and the generators of the 
stabilizer are measured with outcomes a, = tr(pKi), i = 
1, . . . , n. As convention for the coloring we use \A\ Amber 
and \B\ Blue qubits, taking \A\ > \B\. A generator Ki 
is said to be Amber (Blue) if i corresponds to an Amber 
(Blue) qubit. 

Given this tomographically incomplete data, one is 
now interested in finding the minimal entanglement (ac- 
cording to a certain entanglement measure) compatible 
with the measurement data. Mathematically, this is a 
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formulated as the semidefinite program |151 I16| : 

Emin = min p {E(p) : tr(pK. t ) = aup > 0}, (1) 

where E(p) is the entanglement quantifier of choice. We 
will consider the following entanglement measures: the 
relative entropy of entanglement is defined as |17j 

Er(p) = min rTeS Eptr[p{log 2 p - log 2 cr), (2) 

where SEP denotes the set of fully separable states. The 
global robustness of entanglement is given by the mini- 
mum amount of an unnormalized state that has to be 
mixed in to the given state to wash out all entanglement 

R(p) = min a {tr(cr) : p + a £ SEP}. (3) 

Let us return to the estimation of entanglement mea- 
sures from stabilizer measurements. The crucial point is 
that the minimization ([!]) does not need to be carried out 
over all states p. Instead, it suffices to minimize the en- 
tanglement quantifier over stabilizer diagonal states only. 
This can be seen in the following way: since the stabi- 
lizer operators are mutually commuting, the measure- 
ment outcomes tr(pKi) are invariant under any rotation 
of the density matrix of the form p — > KjpKj, Kj G S. 
Due to convexity of the entanglement quantifier, it is le- 
gitimate to apply a local symmetrization procedure , col- 
loquially referred to as "twirling", to the state. This 
is performed by averaging over all stabilizer rotations 
p — > ^a^Zj—\KjpKj. In so doing, the optimization is 
restricted to stabilizer diagonal states of the form 

i 1 

P= ¥l E ' v i ■■■'<»■ (4) 

ii,...,i n =0 

where some coefficients are determined by the measure- 
ment outcomes, while the rest are variables. Since the 
stabilizer operators are mutually commuting, and their 
spectrum is given by {— 1,+1}, it is straightforward to 
compute the eigenvalues of p as Aj = ?c ?- 
Note that a stabilizer diagonal state corresponds to a 
mixture of graph states generated by these stabilizer op- 
erators. This can easily be seen in the following way: 

p = E h\Gd <G t i = E E(- 1 ) ?fc ^ 1 

k k i 

(5) 

= e with c< = i- E(-ir^. (6) 

i k 

Estimating the robustness of entanglement In order 
to estimate the global robustness of entanglement, we 
begin by bounding it from below in the following way: 



for any mixed state p = XkPk and any index m it 
holds that [20] : 

R{p) > M 1 + R (Pm)) ~ 1- (7) 

Since we minimize the robustness over twirled states 
of the form p = ^2%\%\G%){G%\, a lower bound on 
the robustness of entanglement consistent with the sta- 
bilizer measurements is provided by the minimum fi- 
delity F that may be inferred from such measure- 
ments. The minimization of the least fidelity com- 
patible with stabilizer measurements reads: F — 
min p [tr(p\G {kl .... Mn) ){G [ku ... Mn) \) : tr(pK t ) = di,p> 0]. 
By Lagrange duality, one finds the dual problem: F = 

max^iJ^iP-iai : \G(k 1 ,...,k n ))(G(k 1 ,...,k n )\ ~ > °]- 

Solutions to primal and dual can be constructed for an ar- 
bitrary number of qubits, attaining the following analyt- 
ical optimal solution [TB]: F = max[0, a-i — n + 2)]. 
Combining this with Eq. Q provides us with the follow- 
ing lower bound on the global robustness of entanglement 
that can be achieved from stabilizer measurements: 

Rmin(p) >max{0,^ B ^max[0, 1(^a(-n+2)]-l}. (8) 

i 

Estimating the relative entropy In a similar fashion, 
we now calculate the lower bound on the relative entropy 
on entanglement in the case of stabilizer measurements. 
First, note that once more the optimization may be re- 
stricted to stabilizer diagonal states resp. mixtures of 
graph states. Then, the lower bound on the relative en- 
tropy is given by 

ER m i n (p) > max{0, \B\-max[S(p) : tr(pK l ) = ai,p> 0]}. 

(9) 

We can prove this in the following way: first, note that 
the two-coloring divides the system in the two partitions 
A and B. Now one uses the fact that the relative entropy 
is lower bounded by the difference between the entropy 
of system A resp. B and the entropy of the total system 
E3|: 

Er(pab) > max(S(p A ), S(p B )) ~ S(p AB )- (10) 

In our case, we consider only mixtures of two-colorable 
graph states, so that tracing out system A results in a 
maximally mixed state with entropy \B\. Hence, the 
minimization of the relative entropy involves an entropy 
maximization. This can be achieved as outlined in [21] : 
measuring the generators of the stabilizer group gives 
rise to probability distribution pjf^ = 1± 2 ak for the pro- 
jections upon the stabilizer eigenspaces. Furthermore, 
we denote the probability distribution of the joint state 
of the system by A^...^. A crucial feature of the en- 
tropy is subadditivity: S(X il ... in ) < J2k=i,s=± H (Pk)' 
where H denotes the classical entropy function. A 
little thought shows that the above inequality holds 
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with equality for the probability distribution given by 
Aj x ...i„ = Ilfc=i 1+ ^~2~ > ~ ; thus giving the exact max- 
imal entropy S max = ~Y!,i 1 ...i n =a*'h.~iJogXi 1 ...i„- To 
conclude, the lower bound on the relative entropy of en- 
tanglement that can be inferred from stabilizer measure- 
ments is computed as 

1 

ERmin = max{0, \B\ + A, V/A,,..., \. (11) 

i\ ...i n — 

Upper bounds In the previous paragraphs we have de- 
rived lower bounds on the minimal entanglement that is 
consistent with the statistics obtained from measuring 
individual stabilizer operators. It is also of use to derive 
upper bounds. A simple approach to doing this is avail- 
able using the results of [50]. The total Hilbert space 
can be divided up into subspaces, where each subspace 
is labelled by a deterministic outcome for all the Am- 
ber stabilizers. For any pure two colourable graph state 
the entanglement Er is given by \B\, and the robust- 
ness of entanglement is given by 2' B I — 1. In 2D] it was 
shown that for a mixed (twirled) state supported entirely 
in such a subspace, Er is given by Eji(p) = \B\ — S(p), 
whereas R{p) = 2' B I max A^— 1. Let us use the symbol a 
to denote a possible set of outcomes for the Amber mea- 
surements, and let b denote a possible set of outcomes for 
the Blue measurements. Hence any state that is diago- 
nal in the graph state basis can be described as a prob- 
ability distribution p(a, b) = p(a)p(b\a) corresponding to 
the probabilities for getting the various possible stabilizer 
outcomes. We can partition such a state into a mixture of 
states that are individually supported on each of the Am- 
ber subspaces, such that p = ^2p{a)p a , where a is a bit 
string corresponding to the positive/negative stabilizer 
subspaces of the Amber qubits. By concavity of the en- 
tropy function we find that: Er(p) < \B\ — ^2 a p(a)S(p a ). 
But S(p a ) is given by a classical entropy H (p(b\a)) where 
p(b\a) is the conditional probability distribution for get- 
ting outcomes b upon finding a. Thus we obtain: 

E R (p)<\B\-J2p(a)H(p(b\a)). (12) 

a 

Similarly, since the A subspace entanglement is given by 

R(p a ) = 2 lB KKna X (Pa)-l), we find 

R(p) + 1 < 2 |s| ^p(a)max b p(b\a). (13) 

a 

Hence to get upper bounds to the minimal entan- 
glement consistent with the measurement outcomes, we 
need to pick the p(a, b) consistent with the marginal dis- 
tributions that minimises these expressions. The relative 
entropy can now be estimated by noticing that the condi- 
tional entropy is upper bounded by H{p{b)) and choosing 
a product distribution p(b) = p{bi)...p{b\B\), f° r which it 
is well known that it maximises H{p{b)). Thus, we obtain 

E Rmin (p)<\B\-H(p(b)). (14) 



Let b* be the maximiser of p(b) . As b* is a specific value of 
b, it holds that p(a) m&Xbp(b\a) > p(a)p{b) = p(a, b*). It 
follows that J2 a P( a ) maxf,p(6|a) > p(b*), a lower bound 
that is tight as it is attained by setting p(a, b) — p{a)p{b). 
Since we do not have complete information about p{b*), 
we need to minimise it subject to the individual stabilizer 
statistics p(b\),...,p{b\B\)- This problem is equivalent to 
minimising the Zoo norm of a joint probability distribu- 
tion p(b) = p(bi, b\B\) constrained to fixed marginals 
p(bi),...,p(b\B\), which is equivalent to the fidelity min- 
imisation (only considering the Blue stabilizers). Hence 
we obtain: 



R mm {p) + 1 < 2^max 




(15) 

Quality and scaling of the bounds In order to check 
the quality of the entanglement estimates, we consider 
noisy two-colorable graph states. Assuming the exper- 
iments starts from a perfect graph state, which is then 
subjected to local dephasing for a certain time, we can 
take the density matrix time evolution to be governed by 
the following master equation: 

P= -^C^ZipZi- p), (16) 

i 

where 7 is the dephasing constant. The effect of such 
noise on graph states has been studied in detail in Ref. 
|22) . Due to the dephasing the stabilizer coefficients suf- 
fer a decay exponential in the dephasing constant. For 
our test we consider a linear chain of qubits subject to 
this noise. It can be shown that the stabilizer coefficients 
obey the following time evolution in this noise model: 
Ci 1 ...i n (t) = exp{—^t^2 k if t ). Estimates according to the 
described methods for the logarithm of the global robust- 
ness of entanglement and the relative entropy of entangle- 
ment are shown in Figures [l] and resp. [2j The robustness 
of entanglement can be estimated up to a certain num- 
ber of qubits, for which a non-zero fidelity can be inferred 
with the target state. This effectively sets a threshold to 
the estimation. In contrast, the relative entropy may be 
estimated even for larger noisy systems without suffer- 
ing from the threshold problem. However, the difference 
between lower and upper bounds apparently grows with 
system size. 

In many cases, stabilizer measurements are carried out 
via local measurements. This local information on the 
quantum state could in principle be used to improve the 
bounds on minimization of entanglement measures from 
incomplete information on the density matrix, since they 
restrict the set of separable states involved in the opti- 
mization. Note however, that local measurement opera- 
tors generally do not commute with the stabilizer opera- 
tors. This implies that symmetries cannot be exploited. 
For this reason we do not consider local measurement 
data in our scheme. 
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tion from incomplete measurement data usually rely on 
numerical methods such as convex optimization which is 
limited to systems no larger than 20 qubits - even if sym- 
metries can be exploited. Our scheme should therefore 
be invaluable for future graph-state experiments. 

Our results may also be interesting to study the effects 
of various noise models on the entanglement dynamics of 
two-colorable graph states. One step in this direction has 
been made in Ref. [24] , where the entanglement of graph 
states under the influence of Pauli maps is investigated. 

This work was supported by the EU Integrated Project 
QAP and EU STREP projects HIP and CORNER. MBP 
acknowledges an Alexander von Humboldt Professorship. 



FIG. 1: Upper and lower bounds on the estimate of the log- 
arithmic global robustness of entanglement for linear graph 
states subject to local dephasing. A non-zero estimate of the 
entanglement is possible as long as a non-zero fidelity with the 
graph state may be inferred from the stabilizer measurements. 



lower boundy 


= 0.01 


lower boundy 


= 0.03 


lower boundy 


= 0.10 


upper boundy 


= 0.01 


upper boundy 


= 0.03 


upper boundy 


= 0.10 




number of qubits 



FIG. 2: Upper and lower bound on the estimate of the rela- 
tive entropy of entanglement for linear graph states up to 1000 
qubits subject to local dephasing. In contrast to the robust- 
ness of entanglement there is no limit to the lower bound, but 
the difference between upper and lower bounds grows with 
system size. 



Conclusion Here we have shown how entanglement of 
arbitrarily large graph states can be estimated from sim- 
ple measurements. High-quality bounds on the robust- 
ness of entanglement and the relative entropy of entan- 
glement have been derived for stabilizer measurements. 
The stabilizers of two-colorable graph states can be mea- 
sured in two measurement settings (assuming the mea- 
surements can be performed simultaneously), thus our 
scheme avoids the exponential overhead required by full- 
state tomography. In addition, the results presented here 
are of an analytical form that allows for extremely effi- 
cient post-processing. In contrast, quantum state tomog- 
raphy requires computationally hard post-processing of 
the measurement data to create an estimate of the real 
density matrix, and schemes for entanglement estima- 
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